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Abstract. - We calculate the roughness correction to the Casimir effect in the parallel plates 
geometry, for metallic plates described by the plasma model. The calculation is perturbative in 
the roughness amplitude, with arbitrary values for the plasma wavelength, the plate separation 
and the roughness correlation length. The correction is found to be always larger than the 
result obtained in the Proximity Force Approximation. 



High precision measurements of the Casimir force [1] have been reported during the last 
years [2]. They allow for an accurate theory/experiment comparison [3] and, consequently, 
for a search for new weak forces with submillimetric ranges [4] . Theoretical predictions have 
to take into account the differences between realistic models of the experiments and the ideal 
configuration initially considered by Casimir. Most experiments are performed with a plane- 
O" 1, sphere geometry rather than a plane-plane geometry. Temperature corrections to be added 
to the vacuum contribution play an important role when the distance L between the plates is 
above 1/um. Finite conductivity and roughness of the metallic plates used in the experiments 
provide the major corrections for the distances of the order of a few hundred nanometers 
probed by the most accurate experiments. The spatial variations of the surface potential also 
affect the force measurement [5] . 

All these effects must be considered simultaneously since they affect a single observable, 
the Casimir energy. For example, the thermal and plasma corrections cannot simply be 
multiplied in the intermediate range where both effects are noticeable [6] because they are 
in fact correlated to each other. In this letter, we study the correlation arising between the 
conductivity and roughness corrections at the short distances where both are appreciable. 
To this aim, we describe the optical response of the metallic plates by a plasma model with 
a dielectric function e = 1 — ujp/aj 2 where u>p is the plasma frequency. We evaluate the 
roughness effect perturbatively in the roughness amplitude. This effect then depends on the 
hierarchy between the other relevant length scales, namely the plate separation L, the plasma 
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wavelength Ap = 2ttc/lup and the correlation length lc which characterizes the roughness 
spectrum. 

Wc first consider a plane-plane geometry and define the surface profiles by the functions 
hi(x,y) (i = 1,2) giving the local heights with respect to the mean separation L along the z 
direction. These functions are defined so that they have zero averages. We consider the case 
of stochastic roughness characterized by spectra 

<r«(k) =yd 2 re~ ,;k - r (Mr)MO)), » = 1,2. (1) 

The surface A of the plates is supposed to contain many correlation areas, allowing us to take 
ensemble or surface averages interchangeably. The two plates are considered to be made of 
the same metal and the crossed correlation between their profiles is neglected. 

The variation of the Casimir energy Epp is calculated to second order in the perturbations 
hi, leading to the following expression for the roughness correction [7] 

SE PP = J 0G(k),(k) (2) 

Cr(k) = CTn(k) + <7 22 (k). 

With our assumptions, the spectrum <r(k) fully characterizes the roughness of the two plates. 
The correlation length £c is defined as the inverse of its width. The response function G (k) 
then describes the spectral sensitivity to roughness of the Casimir effect. Symmetry requires 
that it only depends on k — |k| . The dependence of G on k reflects that not only the roughness 
amplitude but also its spectrum plays a role in diffraction on rough surfaces [8]. It is only 
at the limit of smooth surface profiles k — > that the effect of roughness may be calculated 
from the Proximity Force Approximation (PFA) [9] by averaging the 'local' distances over the 
surface of the plates. 

In previous discussions of the roughness corrections [7], the sensitivity function G (k) was 
analyzed only in the two cases of short (L <C Ap) and long distances (L > Ap). For the short 
range limit, it was deduced from earlier calculations of Maradudin and Mazur [10]. In the long 
range limit, it was derived from the evaluation by Emig et al. [11] of the effect of corrugation 
of a perfectly reflecting plate. In the present letter, we give the results of a new evaluation of 
G (fc) valid for arbitrary separations L. This evaluation relies on calculations of non-specular 
reflection coefficients associated with rough plates, taking into account the roughness-induced 
coupling between Transverse Electric (TE) and Transverse Magnetic (TM) polarizations. The 
full calculations will be presented in a longer paper [12]. Here, we discuss their results which 
allow us to obtain the roughness correction at any distance, in particular at the intermediate 
distances corresponding to most experiments. 

Before entering this discussion, let us emphasize that all these results may be applied to the 
analysis of the plane-sphere geometry employed in most experiments. To this aim, we use the 
PFA to relate this plane-sphere geometry to the plane-plane configuration, which was taken 
as the benchmark for our perturbative calculation. We thus obtain the relative correction of 
the force Fps in the plane-sphere geometry from that evaluated in (0) for the plane-plane 
geometry [7] 

A ^SJZl = 5 _EfL. (3) 
Fps Epp 

Note that applying the PFA to the study of the sphere-plane geometry only requires the sphere 
radius R to be sufficiently large. Besides the usual requirement R ^> L, it is also necessary 



P. A. MAIA NETO et al: ROUGHNESS CORRECTION TO THE CASIMIR FORCE 



3 




0.2 0.3 
k [nrn -1 ] 



0.4 



Fig. 1 - Variation of G/Epp versus k for the distances L — 50nm (solid line), L — lOOnm (dashed- 
dotted line), L = 200nm (dotted line), and L = 400nm (dashed line) for Xp — 136nm. 



to assume RL 3> so that many correlation areas are included in a given nearly-plane 
local section of the sphere. In contrast, using the PFA to calculate the roughness correction 
requires the correlation length to be larger than the separation Iq ^> L. This is clearly a 
more restrictive condition. In the present work, we assume the PFA validity conditions to be 
obeyed for plane-sphere geometry but not necessarily for roughness. 

According to @ , the relative roughness correction is obtained by integrating the ratio 
G (k) /Epp over the roughness spectrum a (k). This ratio is plotted on Fig. 1 as a function of 
the roughness wavevector k for several different values of the distance L. As for all numerical 
examples considered below, we take Xp = 136nm which corresponds to gold covered plates. 
As expected, the relative roughness correction is larger for shorter distances. In the following, 
we discuss the values of G in the limit k — > corresponding to the PFA. We then come to the 
main result of this work, that is to say the ^-dependence of G which reveals the departure of 
the effect of roughness from its PFA description. 

The PFA result is recovered as a consequence of the following identity 



G(k^0) = 



E'pp(L) 



(4) 



where the derivative is taken with respect to the plate separation L. This identity is obeyed 
by our result for arbitrary values of L and Xp [13]. If we now suppose that the roughness 
spectrum <r(k) is included inside the PFA sector where G (k) ~ G (0), G may be replaced by 
its limiting value J1J and factored out of the integral J2J) thus leading to the PFA expression [7] 



A= E'ML) a2 



d^k 

4^2 



(5) 



ff (k) = (hi 



hi). 



In this PFA limit, the correction depends only on the variance a 2 of the roughness profiles, 
that is also the integral of the roughness spectrum. 
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Fig. 2 - Variation of p versus k for several values of L (same conventions as on Fig. 1). 



In the general case in contrast, the sensitivity to roughness depends on the wavevector k. 
This key point is emphasized by introducing a new function p which measures the deviation 
from the PFA [7] 

P{k) = Wy (6) 

This function is plotted on Fig. 2 as a function of k for several values of L. It is almost 
everywhere larger than unity, which means that the PFA systematically underestimates the 
roughness correction. The inlet shows p for small values of k where the PFA is a good 
approximation for the shortest distances, for example L = 50nm. To give a number illustrating 
the deviation from the PFA, we may say that p ~ 1.6 for L = 200nm and k = 0.02nm _1 , 
which means that the exact correction is 60% larger than the PFA result for this intermediate 
separation and a typical roughness wavelength 27r/fc ~ 300nm. 

Fig. 2 indicates that p(k) grows linearly for large values of k. This is in fact a general 
prediction of our full calculations for arbitrary values of L and A p : 

p(k) = ak for fc> LU P /c,l/L. (7) 

The dimensionless parameter a/ L depends on Kp — ujpL/c = 2irL/\p only, and is given by 

HA 



a = 



(2^i> G(0) J Q 
f 2(K 2 



dKK 



K 



K 2 

2n 2 + k 2 



(8) 



x < Kf TE 



n 



2\2 



k 2 {2k 2 - zn 2 



(KK t ) 2 - (k 2 - n 2 



TM 



The dimensionless integration variable K represents the z component of the imaginary wavevec- 
tor multiplied by L [14], and K t — y/ K 2 + Kp corresponds to the z component of the imagi- 
nary wavevector inside the metallic medium. We denote similarly Q the imaginary frequency 
multiplied by L/c. f e are the loop functions describing the optical response of the cavity for 
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the two orthogonal polarizations e = TE, TM: 



= r 2 exp(-2K) 
Je 1 - r 2 exp(-2i?r) ' ( ' 



The corresponding reflection coefficients are given by 



^ = -^1 do) 



( 



K t + K 



K + K, 



In Fig. 3, we plot the coefficient a as a function of L, still with the plasma wavelength of 
gold Xp — 136nm. At the limit of short distances, we recover from (JHJ) our previous result [7] 

a = 0.4492L for fc _1 «L« X P . (11) 

At the limit of large distances, the angular coefficient saturates, yielding 

a = ^-\ P for £T 1 <A P <L. (12) 

o(j7T 

This result is derived from by expanding the integrand in its righthand side in powers of 
Xp around Xp — 0. Remarkably, it differs from the long distance behavior reported in Ref. [7], 
which was derived from the analysis of corrugation for perfectly reflecting plates [11]. In fact, 
the high-fc expression © holds when the roughness length scale 1/fc is much smaller than 
both Xp and L. A different result is obtained when Xp rather than 1/fc is the shortest length 
scale. In order to see it, we derive from our general result an expression for G valid at this 
limit: 

HA 1 f°° dKe~ 2K f K+q , , . 

G{k) =-7-,TW-J 7— / dK' (13) 



{KK') 2 + \{K 2 + K' 2 ~q 2 ) 2 



1 



for X P -> 0. 



Numerical integration of l|13|) agrees with the expression of Emig et al. [11] for arbitrary 
values of q = kL; K has the same meaning already discussed in connection with (JSJ while K' 
corresponds to the longitudinal component of the imaginary wavevector associated with the 
diffracted wave. 

In order to discuss the regime Xp <^ 1/fc <^ L, we now take the high-fc limit of the right- 
hand side of i|13|) . Due to the presence of the exponential factor exp(— 2K), the dominant 
contribution comes from the corner K ~ 1, K' ~ q of the rectangle associated to the integration 
region. We may thus neglect exp(— 2K') and recover the long distance limit of [7]: 

2 HAq f°° K 3 e~ 2K tt 2 a q 

p = -Lk for A P <fc _1 <i. (14) 

In summary, the long-distance behavior is given by H12|) when 1/fc <C Ap < L, and by l|14l) 
when Ap 1/fc <C L. The cross-over between these two regimes is shown in the inlet of 
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Fig. 3 - Variation of the angular coefficient a versus L for Ap = f 36nm. The analytical result for 
k" 1 <C L <C Ap is shown as the dotted line and for k^ 1 <C Ap <C L as the dashed line. A comparison 
between this second result (dashed straight line) and the exact p(k) (solid line) is shown in the inlet 
for L — 2^tm. The analytical result p = Lk/3 predicted by the model of perfect reflectors (dotted 
line) is valid only in the intermediate range Ap < fc" 1 <L. 

Fig. 3, where we plot p as a function of k for L — 2/im. The model with perfect reflectors fails 
when 1/k <C Ap because Fourier components with K ~ 1 <C K p , for which the plates behave 
as perfect reflectors, are diffracted into components with K' ~ q 3> K p which are no longer 
perfectly reflected by the plates. 

In order to reach definite conclusions about the accurate evaluation of the roughness cor- 
rection, the first crucial step is to measure the roughness spectrum cr(k). Once in possession 
of this experimental input, one may calculate the roughness correction from the second-order 
perturbation formula PJl. In order to go further in the present letter, we consider the simple 
model of a Gaussian roughness spectrum [10] 



a 2 is the roughness variance and tc the correlation length. Using this model, we now illustrate 
the preceding results by giving scaling laws obtained in some limiting cases. 

As we increase the distance L between two given plates, characterized by the length scales 
Ap and £c, the roughness correction decreases according to power laws that may be derived 
from Eqs. J2J and l|15|l. We first consider the case of very smooth surfaces Ap -C tc- The 
correction at short distances L -C Ap <C tc may thus be calculated from the PFA Eq. (0) and 
leads to A = 3a 2 /L 2 . The PFA still holds when L is increased beyond the plasma wavelength 
into the intermediate range Ap <C L <C £c giving A = 6a 2 /L 2 . As we increase the distance 
further, the correction decreases at a slower rate. When Ap <C £c "C L, we find from J5J, 
(II4II and 1)15(1 the power law A = 2y / 7r a 2 /(£cL), which represents a correction larger than 
the PFA result by a factor Ljic I. The last two results may be derived from the study of 
perfect reflectors [7,11]. 




(15) 
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On the other hand, if we start with a pair of very rough surfaces tc <C Ap, we find a 
completely different behavior as the distance L is increased [15]. The short distance limit is still 
governed by the PFA power law 1/L 2 but in the intermediate range Iq <C L <C Ap we now find 
from Eq. lfTT|) A = 2.7^/tt a 2 j{lcL). For very long distances lc <C Ap <C L, saturation leads to 
a faster decrease of the correction, and from Eq. (|12[) we find A = (14/5-v/tt) (Ap/^c) { a2 /L 2 ). 
This is a 1/L 2 decrease, as in the PFA, but with an additional factor of the order of \pjlc 3> 1- 

In conclusion, we have computed the second-order response function G{k) for arbitrary 
values of the plasma wavelength Ap, distance L and roughness wavevector k. This allows for 
a reliable calculation of the roughness correction up to second order in the profiles h\ and hi. 
We have derived analytical results in some limiting cases. In particular, we have discussed 
the limits of long distances and short roughness wavelengths and shown that their relation to 
the model of perfect reflectors is much richer than previously thought [7,11]. The PFA results 
are recovered as the limiting case of very long roughness wavelengths. Our present analysis 
proves that they systematically underestimate the roughness correction [16]. This has to be 
taken into account in the search for bounds on new weak forces in the submillimeter range [4] . 

We thank Cyriaque Genet and Marc-Thierry Jaekel for discussions. PAMN thanks Insti- 
tuto do Milenio de Informagao Quantica and CNPq for partial finantial support. 
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